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NOMENCLATURE 

A  -  constant  in  regression  rate  equation  (7) 

B  -  blowing  parameter 

C  -  specific  heat  at  constant  pressure 
P 

E  -  activation  energy 

h  -  enthalpy 

hf  -  heat  of  formation 

k  -  turbulence  kinetic  energy 

m  -  mass  flow  rate 

M  -  average  molecular  weight 

m.  -  mass  fraction  of  species  i 

p,P  -  pressure 

q"  -  heat  transfer  flux 

r  -  radial  distance 

r  -  fuel  regression  rate 

R  -  universal  gas  constant 

R.  -  reaction  rate  of  species  i 

T  -  temperature 

u  -  axial  velocity 

v  -  radial  velocity 

w  -  tangential  velocity 

r.  -  effective  transport  coefficient  of  species  i  =  y/a. 

5  -  incremental  distance  from  wall 

e  -  turbulence  dissipation  rate 

u  -  effective  viscosity 

p  -  density 

a.  -  Prandtl    or  Schmidt   number  for   species      i 

<J>  -  mixture  fraction    (mass   fraction   of  total    fuel) 
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I.   INTRODUCTION 

The  solid-fuel    ramjet    (SFRJ)   most   often   consists   of  a  solid-fuel    grain 
which  provides   the  walls   for  the  combustion   chamber.     A  sudden  expansion   at 
the  air   inlet    (either  axial    or   side-dump)   can   be  used  to  provide  flame 
stabilization.     Downstream  of  the   flow   reattachment   a  turbulent   boundary   layer 
develops   and   includes   a  diffusion-controlled   flame  between  the   fuel-rich   zone 
near  the  wall    and  the   oxygen-rich   central    core.     Due  to  that   diffusion   flame, 
heat   is   transferred   by   convection   and   radiation   to  the   solid   surface,   causing 
vaporization   of  the   fuel. 

At  the  Naval    Postgraduate  School    both  mathematical   modeling*^  and 
experimental    efforts^-lO  have   been   conducted  to  determine  the  effects   of 
design   and  operational    variables   on   the  obtainable  performance. 

The  computer  model    simulation   of  the  SFRJ   combustion   process   has   evolved 
from  an   original    stream-function   vorticity   formulation*   to  a 
primitive-variable    (pressure,   velocity)  model   which   includes   an  aft  mixing 
chamber. 3     These  axisymmetric  models   have  not   included   radiative  heat  transfer 
to  the   fuel    surface.     This   has    limited  the  utility  of  the  model,   since  many 
all -hydrocarbon   fuels   produce   significant   amounts   of   radiative  transfer 
through  the  generation   of   carbon   particulates    in  the   flame  zone.     Recently'*,   a 
simplified   radiative  model    was    included   in   the  primitive-variable  model    and 
resulted   in   improved  prediction   of  the  fuel    regression-rate  profile.     The 
primitive-variable   codes  were   based  on   the  Champion   2/E/FIX  computer  program 
developed  by   Pun   and  Spalding^.     The  Champion   code   is   a  two-dimensional    code 
and  therefore,   could  not   be  used  to  model    two  important   SFRJ   geometries; 
side-dump  and  bypass.     The   side-dump  configuration   has   air  injected  through 


the  side   of  the   fuel    grain   and  the  dome   region   is   constructed  of   fuel.     The 
bypass   configuration   ducts   a  portion    (typically  50%)   of  the  total    air  flow 
into  an   aft  mixing  chamber.     A   full    three-dimensional    (3-D)   model    is   needed  to 
predict  the  flow-field   in  these   configurations. 

The  side-dump  configuration   is   a  very  difficult   one  to  develop/optimize 
using  only  experimental    data.     There   is   a  complex   interaction   between  the 
location/size   of   the   inlet   dump,    the   length   of  the  dome   region   upstream  of  the 
dump,   and  the   resulting   fuel    regression    rate  pattern.     Regression   rate  data 
for  the  side  dump   geometry  have  not   appeared   in  the  open   literature.     However, 
several    general    observations   can   be  made.      Swirling  the   inlet   flow  can 
increase  and   smooth-out   the  downstream  fuel    regression    rate  pattern. 
Regression   rates  within  the  dome   region,   upstream  of  the  dump,   can   increase  or 
decrease,    depending   upon   the  dome    length,    number   of   inlets   and  amount   of 
swirl.      3-D  models   are,    in  this   case,   a  necessity   in  the  combustor  development 
process. 

Several    3-D,   elliptic,   primitive-variable  codes   have  been  developed.  .  One 
such   code,   developed   by  the  AiResearch  Manufacturing  Company, 12, 13   -jS   an 
extension   of  the  Champion   code.     The   former  computer  code  was   developed  to 
model    gas   turbine   combustors   and   formed   the   basis   of  the   3-D,   SFRJ  model 
discussed  herein. 


II.  MODEL  OVERVIEW 

The  Garrett/AiResearch  model  is  3-D,  elliptic,  for  steady,  subsonic  flow 
and  includes  finite-rate  chemical  kinetics.  The  dependent  variables  are 
u,  v,  w,  h,  k,  e,  <(> ,  n,  ,  m-.,,  ra,  and  mr~.   In  addition,  if  radiation  is 

included  (a  six-flux  model),  then  three  radiation  flux  vectors  are  calculated. 

Soot  emissions  can  be  calculated  if  desired.  Liquid  spray  calculations  can 

also  be  made  but  were  not  necessary  in  the  SFRJ  model.  The  turbulence  model 

is  the  same  as  used  in  the  Champion  program,  namely  the  modified  two-equation 

(k-e)  model  of  Jones,  Launder  and  Spalding^"^. 

Hydrocarbon  combustion  is  considered  to  be  a  four-step  process  as 

fol lows : 

C  H  +  C  H  ,  +  H,  (1) 

x  y    x  y-2    2  v  ' 

C  H  „  +  4  0„  -►  xCO  +^  H0  (2) 

x  y-c       2  2         2   2  v  ' 

CO  +  |  o2  -►  co2  ,2) 

H2  +  j   02  +  H20  (4) 

Arrhenius  rate  expressions  are  used  for  each  of  the  above  reactions.  These 
chemical  kinetics  limited  rate  expressions  were  compared  to  rate  expressions 
based  on  an  eddy-breakup  model.  The  latter  expressions  attempt  to  account  for 
the  effects  of  turbulence  intensity  and  scale  and  species  concentration  on  the 
consumption  rate  of  a  particular  reactant.  The  two-part  boundary  layer 
(divided  by  y+  =  11.5)  used  in  the  Champion  code  is  also  employed.  Details  of 
the  computer  program,  the  differential  equations  and  the  line  by   line,  finite 
difference  solution  procedure  are  presented  in  references  12  and  13. 


III.  MODEL  MODIFICATIONS  AMD  SOLUTION  PROCEDURES 

The  major  modifications  required  of  the  Garrett/AiResearch  model  were 
geometric  and  the  inclusion  of  a  fuel  wall  on  .vhich  finite-rate  chemical 
kinetics  could  occur.  The  geometric  changes  required  to  model  the  SFRJ 
configurations  were  readily  incorporated.  Several  additional  variables  and 
subroutines  were  necessary  to  include  the  regressing  fuel  surface  effects.  In 
addition,  fuel  property  subroutines  were  changed  to  reflect  the  properties  of 
Plexiglas  (PMM)  and  HTPB,  the  two  solid  fuels  used  in  the  present  study. 

The  governing  equations  on  the  boundary^  for  unity  Lewis  number  are: 

Energy:    q  ".  =  (pv)  (h  -hj  +  r.   (—•)  (5) 

*J  rad   v  'wv  w  f         h  v3ry  w  v  ; 

w 

(q  can   be   calculated  within   the  existing  model    but  was 

vVad  3 

neglected  in  the  present  investigation.) 

3m. 

Species:        R.      =r.    (^-1)       +   (pv)    (m.    -m.    )  (6) 

f  i  i    v9  r    'w  v      'wv    l        i    y  K    J 

ww  w       g 

(the   values   of  R-j   were  calculated   in   the  same  manner  as   for  the 

internal    nodes) 

It  was   assumed   that   the   fuel    regression    rate   could  be   represented   in   an 

Arrhenius   format:^ 

Pfr  =    (pv)w  =  Apfe"'     w  (7) 

Using  the  Couette   flow  approximation   for  the  boundary   layer  behavior  with  mass 
transfer*'^ 

u  *n(l+B)  (8) 

1  B=0       B  l    ; 


where  B  can  be  evaluated  from  the  enthalpy  as 

h  -h 
_P 

R  -h.  -  (q  "./(pv)  ) 
w  fg  VMrad/x   ;wy 


B  -  t: — c nr- n — tttt,  \ — r  (9) 


The  enthalpy  can  be  expressed  as 

T 

P 


h   =  E  m, (h"  +  J  W  C  dT)  (10) 


ref 

The  solution  procedure  for  the  boundary  conditions  on  the  fuel  wall  was 
as  follows: 

(a)  estimate  T  and  B 

\  /  w 

(b)  calculate  (pv)   from  (7) 


(c)  calculate  r.  and  r   from  (8) 

w      w 

(d)  calculate  n.   for  all  species  from  (6) 


(e)  calculate  the  other  species  (mH  ,,,  mnx,  mrn  ,  r\.  )  from 
species  conservation  equations      2         2    2 

(f)  calculate  hw  from  (10) 

(g)  calculate  B  from  (9)  with  q  "  =  0 
Iterate  (c)  to  (g)  until  B  converges 

(h)  calculate  (pv)   from  (5) 

(i )  calculate  Tw  from  (7) 

Iterate  (b)  to  (h)  until   (pv)w  converges 

( j )  calculate  pw  from 

pw  =fl 
H 

(k)  calculate  v   from   ^pv^w 
v  '  w       

pw 


It  was  found,  as  previously  found  by  Netzer*  and  Stevenson  and  Pletzer^, 
that  the  near-wall  turbulence  dissipation  had  to  be  increased  on  the  step  face 

(e  =k  3/2/0.46)  and  that  the  grid  spacing  adjacent  to  the  fuel  surface  had  to 

be  maintained  slightly  less  than  that  required  for  y+  =11.5  in  order  to 

P 
obtain  temperature  distributions  in  qualitative  agrement  with  experimental 

data. 

In  order  to  obtain  convergence  it  was  necessary  to  input  initial 
conditions  which  were  in  qualitative  agreement  with  the  final  solution  (for 
example,  the  fuel  mass  fraction  had  to  be  specified  as  decreasing  from  the 
wall)  and  under  relaxation  (typically  0.1  to  0.5)  had  to  be  used. 

The  accuracy  of  the  solution  for  the  entire  combustor  of  the  solid  fuel 
ramjet  is  particularly  sensitive  to  the  boundary  conditions  on  the  fuel  wall. 


IV.  RESULTS  AND  DISCUSSION 

Four  geometries  were  investigated  (Fig.  1);  axisymmetri c  with  axial 
inlet,  bypass  with  axial  upstream  inlet  ,  single  side-dump  and  dual  side-dump. 
All  cases  included  an  aft  mixing  chamber.  For  the  axi symmetric,  axial-inlet 
geometry  the  predicted  reattachment  "points"  (where  u  near  the  wall  was 
interpolated  to  zero)  are  shown  in  Fig.  2.  The  reattachment  "points"  were  in 
good  agreement  with  the  maximum  heat  transfer  locations  measured  by  Krai  1  and 
Sparrow^  but  shorter  than  those  measured  by  Phaneuf  and  Netzer°  in 
nonreacting  flow  with  wall  mass  addition. 

Normalized  regression  rate  profiles  are  shown  in  Fig.  3.  The 
experimentally  measured  profiles  for  PMM  fuel  were  very  sensitive  to  inlet 
turbulence/distortion.  However,  the  3-D  code  (with  finite  rate  kinetics) 
predicted  results  very  nearly  identical  to  the  mixing  limited,  2-D  code^, 
i.e.,  the  maximum  regression  rate  was  predicted  to  occur  upstream  of  the 
experimentally  measured  point.  The  average  regression  rate  was  in  good 
agreement  with  experiment. 

Examples  of  typical  radial  profiles  for  u,  mn  ,  mN  ,  and  T 

near  the  grain  exit  are  shown  in  Fig.  4.  The  effects  of  including  finite-rate 
kinetics  and/or  the  injection  of  O2  that  occurs  into  the  boundary  layer  near 
the  reattachmnt  zone  are  readily  apparent;  a  broad  zone  of  high  temperature 
near  the  maximum  value  and  finite  concentrations  of  oxygen  below  the  flame. 

Figure  5  shows  regression  rate  profiles  at  three  axial  locations  for  the 
single  side-dump  geometry  without  inlet  air  swirl.  Non-uniform,  high 


regression  rates  occur  in  the  dump  region.  Lower  regression  rates  in  the  dome 
region  and  more  uniform  regression  rates  near  the  grain  exit  are  also 
observed. 

Figure  6  shows  the  predicted  regression  rate  results  with  two,  180° 
opposed  side-dumps,  with  and  without  inlet  air  swirl.  The  effects  of  inlet 
swirl  on  this  particular  geometry  are  readily  apparent.  Inlet  swirl  decreased 
the  regression  rates  in  the  dome,  but  significantly  increased  the  rates 
downstream.  Regression  rate  uniformity  was  also  significantly  improved  except 
adjacent  to  the  inlet  dump. 

Figure  7  shows  the  results  of  increasing  the  dome  length  and/or  the  inlet 
dump  area.  Increasing  the  inlet  area  resulted  in  increased  fuel  consumption 
in  the  dome  region  and  increased  nonuni formity  downstream.  Lengthening  the 
dome  above  the  reference  value  significantly  reduced  head-end  fuel  regression 
rates  without  major  effects  in  the  downstream  regions. 


V.  CONCLUSIONS 

The  model  predictions  show  that  significant  variatiions  in  fuel 
regression  rate  (both  ci rcumferential ly  and  axial ly)  are  caused  by  the  manner 
in  which  the  air  is  injected  into  the  fuel  grain.  In  general,  inlet  swirl 
increases  and  smooths-out  the  fuel  regression  rate  downstream  of  a  side-dump 
inlet.  The  effects  of  inlet  location  and  size  on  the  regression  rate  behavior 
are  difficult  to  estimate  a-priori  for  a  specific  set  of  test  conditions.  The 
model  must  be  used  to  predict  the  expected  results.  Much  additional  work  is 
required  in  order  to  validate  the  predictions  of  the  3-D  SFRJ  model.  However, 
the  ability  to  examine  the  effects  of  inlet  air  location  and  flow 
characteristics  on  the  fuel  regression  rate  pattern  (and  the  combustion 
efficiency,  etc.)  has  provided  a  needed  improvement  to  the  earlier  2-0  model. 

Current  work  is  being  directed  at  incorporation  of  the  soot 
production/consumption  and  radiation  subroutines  into  the  SFRJ  configuration 
and  combustion  environment.  Regression  rate  profiles  for  the  more  recent 
configurations  are  also  being  compared  to  experimental  data  as  they  become 
avai 1  able. 
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